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Abstract 

CN ' In this paper we study the asymptotic phase space energy distribution of solution of the 

Schrodinger equation with a time-dependent random potential. The random potential is 
assumed to be with slowly decaying correlations. We show that the Wigner transform of a 
solution of the random Schrodinger equation converges in probability to the solution of a 
radiative transfer equation. Moreover, we show that this radiative transfer equation with 
long-range coupling has a regularizing effect on its solutions. Finally, we give an approxi- 
mation of this equation in term of a fractional Laplacian. The derivations of these results 
' are based on an asymptotic analysis using perturbed-test-functions, martingale techniques, 

\ and probabilistic representations. 
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Introduction. 

t— i ; 

The Schrodinger equation with a time-dependent random potential has been studied for a long 
£T) • time and arises in many application [IJJ HTJ Q21 121 EH ESj , for instance in wave propagation 

in random media under the paraxial or parabolic approximation El El El [TU] . The present 
work has been motivated by data collections in wave propagation experiments showing the 
possibility to encounter propagation medium with slowly decaying autocorrelation functions 
[TBI EE] • These observations have stimulated this field of research [221 EH E21 E7J , but mainly 
in one-dimensional propagation media. These media are very convenient for mathematical 
^ ■ studies, but not relevant in many applications. Random media with long-range correlations 

have been also considered in motion of particles. In this contexts, the authors have shown that 
the deviation of the trajectory of a particle from its mean trajectory converges to a fractional 
Brownian motion [2U1 12B1 127] . 

Let us consider the random Schrodinger equation 

id t 4> + « - VtV(t, x)0 = 0, t > and x G R d , 
0(0, x) =0 o (x), 

with a random potential V(t,x), which is a spatially and temporally homogeneous mean-zero 
random field. Here, t > represents the temporal variable, x £ M. d the spatial variable with 
d > 1, and e <C 1 is a small parameter which represents the relative strength of the random 
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fluctuations. A classical tool to study the phase space energy density of solutions of the random 
Schrodinger equation is the Wigner transform defined by 

W(t, x, k) = J dye ik ^(i, x-f>(t,x + !). 

We refer to |23t 12"5] for the basic properties of the Wigner transfom. In our problem the size 
of the random variations are small, so we have to wait for long propagation distance and large 
time propagation to observe significant effects. Consequently, we consider the rescaled field 

e (t,x) = <£(~,~) 

which satisfies the scaled random Schrodinger equation 

ied t (j) e + —A x <j> e - VeV(-, -)d) e = 0, t > and x e R d , 
2 Ve e / 

0e(O,x) = 0O,e(x). 

In this regime, we consider the scaled Wigner transfom given by 

W e (t,x,k) = J dye ik y^(t,x-e|)<Ae(t,x + e|). 

The scaled Wigner transform is well suited to study functions oscillating on scales of order e _1 , 
since the difference between cf) t (t,x — ey/2) and 4> € {t, x+ey/2) is of order 0(1). In several paper 
[SI El EZ1 Q2J [13 [Ml 133 BE] it has been shown that the expectation of the Wigner transform 
E[W e (t,x, k)] converges as e goes to to the solution W of the radiative transport equation 

d t W + k-V^W = J dpa(p,k)(W(t,x,p) -W(t,x,k)), 

where the transfer coefficient <r(p, k) depends on the power spectrum of the two-point correla- 
tion function of the random potential V. This result holds under mixing assumptions on the 
random potential V. Moreover, it has been shown that the limit W is often self-averaging, that 
is, W e converges in probability to the deterministic limit W for the weak topology on L 2 (M. 2d ) 

In [10] the authors study the Schrodinger equation with a random potential with either 
rapidly or slowly decaying correlations. This paper constitutes a first step in the study of wave 
propagation in random media of dimension strictly greater than 1 with long-range correlations. 
In [TU] the authors study the field <j> itself and not its phase space energy density. In the rapidly 
decorrelating case the authors show that the field and the phase space energy density evolve at 
the same scale, which is of order e~ x . As a result, in this case the scale e _1 is universal in the 
sense that it does not depend on the random potential V . However, in the slowly deccorelating 
case, the authors have observed macroscopic effects happening on the field <j> at a shorter scale 
£ -i/(2«o) j w ith k > 1/2. On this scale, the asymptotic field has a random phase modulation 
given by a fractional Brownian motion with Hurst index kq. Let us note that the parameter 
kq depends on the statistic of the random potential V, and therefore the scale at which we can 
observe significant effects on the field <p is no longer universal. Moreover, as we will see in this 
paper, at the scale e _1 ^ 2K °) the phase space energy of the field </> e is not affected. In fact, the 
phase space energy distribution in both cases, rapidly and slowly decaying correlations, evolve 
on the scale e _1 . 

In this paper, we investigate the propagation of the phase space energy density of the field 
(j) t in a random medium with slowly decaying correlations. We show in Theorem 12.21 that the 
Wigner transform of cf) € converges in probability, for the weak topology on L 2 (M. 2d ), to the 
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unique solution of a radiative transfer equation similar to the one obtained in [B] under rapidly 
decaying correlations, but with an important difference coming from the long-range correlation 
assumption. In contrast with the rapidly decorelating case jnj, the scattering coefficient S(k) = 
J dp<r(k, p) = +00 is not defined anymore. The radiative transfer equation is still well defined 
because of the difference W(t, x, p) — W(t, x, k) which balances the singularity given by the 
long-range correlation assumption. This result shows a qualitative and thorough difference 
between the two cases, rapidly and slowly deccorelating cases. First, in contrast with the 
rapidly decorrelating case, for which the phase and the phase space density evolve on the 
same scale e , the phase of <f> and its phase space energy evolve at different scales in the slowly 
decorrelating case, which are respectively e _1 /( 2K o) and e _1 . Moreover, we show in Theorem 13. II 
that the long-range correlation assumption implies a regularizing effect of the radiative transfer 
equation. This regularizing effect is a consequence of S(k) = +00, and cannot be observed 
under rapidly decaying correlations. Finally, we give in Theorem 15.11 an approximation of 
the radiative transfer equation in term of fractional Laplacian, which permits to exhibit from 
its solutions a damping coefficient obeying a power law with exponent lying in (0,1). This 
approximation corresponds to the long space and time diffusion approximation of the radiative 
transfer equation. 

The organization of this paper is as follows: In Section[U we present the random Schrondinger 
equation that will be studied in this paper. Then, we present the construction of the random 
potential that will be used to model the random perturbations. Finally, we introduce the 
long-range correlation assumption used throughout this paper. In Section 2, we state the main 
result of this paper. We describe the asymptotic evolution in long range random media of the 
phase space energy density of the solution of the random Schrodinger equation. In Section [31 
we present the regularizing effect of the radiative transfer equation. In Section [H we present 
the probabilistic representation of the radiative transfer equation. In Section [5] we present 
the approximation of the radiative transfer equation in term of fractional Laplacian. Finally, 
Sections El El El and [9] are devoted to the proofs of Theorem 12.21 Theorem 13.11 Proposition 14.11 
and Theorem 15.11 

Acknowledgment 
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1 The Random Schrodinger equation 

This section introduces precisely the random Schrodinger equation that we study in this paper. 
We consider the dimensionless form of the Schrodinger equation on M rf with a time-dependent 
random potential: 



with 7 G (0, 1). The case 7 £ (0, 1) permits to address a more convenient study of the phase 
space energy density of the field <ft. The case 7 = 0, which will not be addressed in this paper, 
leads to much more difficult algebra. In ([I]) the strength of the random perturbation is small, 
so we consider the rescaled field 



to observe significant effects after a sufficiently large propagation distance and propagation 
time. Therefore, the rescaled field 4> t satisfies the scaled Schrodinger equation 




(1) 




ied t 4> e + 7rA x e 



1-7 
— e 2 



y (^iT7'7H = with &(0,x)=(Mx) 



(2) 
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with 7 G (0, 1). Here A x is the Laplacian on M d given by A = Y?j=i d% ■ (V(t, x), x G R d , t > 0) 
is the random potential, whose properties are described in the next section. Moreover, the initial 
datum <^o,e(x) = </>o, e (x, C) is a random function with respect to a probability space (S, fj,(dC)), 
and independent to the random potential V. This randomness on the initial data is called 
mixture of states. This terminology comes from the quantum mechanics, and the reason for 
introducing this additional randomness will be explained more precisely in Section 11.31 



1.1 Random potential 

In this section we present the construction of the random potential V which is considered 
throughout this paper. It is also a short remainder about some properties of Gaussian random 
fields that we use in the proof of Theorem 12.21 All the results exposed in this section can be 
shown using the standard properties of Gaussian random fields presented in (TJ [2] for instance. 

In this paper, the random perturbations are modeled using a stationary continuous random 
process in space and time, denoted by (V(t, x),i > 0,x G M d ). We construct our potential 
in the Fourier space as follows. Let Ro be a nonnegative function with support included in a 
compact subset of M. d containing 0, such that .Ro G L 1 (M d ), Rq(— p) = Ro(p), and Ro has a 
singularity in 0. Let us consider 

H = \ ip such that / dpRo(p)\ip(p)\ 2 < +°° r > 

I JR d J 

which is a Hilbert space equipped with the inner product 

(ip,^) H = J dp r (pMp)W) v(y?,vO e n 2 . 

Let us consider (V(t, -))t>0 be a stationary continuous zero-mean Gaussian field on %' with 
autocorrelation function given by 

E[V(h,dpi)V(t 2 , dp 2 )] = {2Ti) d R{t 1 - t 2 , PiWPi + P2) 

and 

E[V(t 1 ,dp 1 )V(t 2 ,dp 2 )} = (27r) d R(t 1 - t 2 ,px)5( Pl - p 2 ), 

and where %' is the dual space of H. In other words, Vn G N*, V(</?i, . . . , <p n ) G % n and 
V(ti...,t n ) G [0,+cx)) n , 

{(V(tl),<pl) W;H (V(t n ),<Pn) H r )H ) 

is a zero- mean Gaussian vector with covariance matrix given by: V(j, /) G {1, . . . , n} 2 



E 

and 

E 



(V(t j ),ip j ) H , jH (V(t l ),ip l ) n ,^ = J^dp p j (p)( Pl (-p)R(t 1 -t 2 ,p) 



(y{ij) t ^j) WH (y{ti),ipi) WH = /_ dp ip j (p)ipi(p)R(t 1 -t 2 ,p). 



Here, the spatial power spectrum is given by 

R(t,p) = e-^R (p), (3) 

where the nonnegative function q is the spectral gap. We assume that the spectral gap is 
symmetric, that is g(p) = q(— p). Particular assumptions involving the spectral gap q will 
be introduced at the end of this section to ensure the long-range correlation property of the 
potential V. 

According to the shape of the autocorrelation function ([3D, we have the following proposi- 
tion. 
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Proposition 1.1 Let 

Tt = <r(V(s,-),s<t) (4) 
be the a -algebra generated by (V(s, •), s < t). We have 

E[V(t + h,-)\T t ]=e-^ h V(t,-) (5) 

andV(ip,ip) G U 2 



E 



(V(t + h), <f) Wjn (V(t + h),il>) Wfll -E[(V{t + h), <p) n , tH \?tM(V(t + h),i/i) WtH \Ft] 

= Jdp ^(p)V(-p)^o(p) (l - e- 2 ^ h ) . 



(6) 



These two properties will be used in the proof of Theorem l2.2[ which is based on the perturbed- 
test-function method. 

Let us note that {V,^)^,^ is a real-valued gaussian process once <p E H satisfies </?(p) = 
(p(— p). According to this last remark, let us introduce the real random potential V defined by 

F(i,x) = (V(t,-),e*) n , jn = j±y d J R V(t,d P )e^, 

where e x G H is defined by e x (p) = e* px /(27r) d . Consequently, the random potential V is a 
stationary real- valued zero-mean Gaussian field with a covariance function given by: V(iij £2) G 
[0, +00) 2 and V(xi, x 2 ) G R 2d 

J R(ti-t 2 ,x 1 -x 2 )=E[y(t 1 ,x 1 )y(t 2 ,x 2 )] / dpii(ti-t2,pK p - (xi - X2) 

= / ^dp J R(w,p)e icj(tl -* 2) e ip -( Xl - X2 \ 

(7) 

where 

fl(c,p) = ^|M ( 8 ) 
w 2 + g 2 ( P ) 

According to the previous construction and [2 Theorem 2.2.1] the random potential V is 
continuous and bounded with probability one on each compact subset K of R x M d . This fact 
comes from the relation 

E[(V(t, x) - V(s, y)f] 1/2 <c(J dp^o(p)) (|t - s\ + |x - y|) , 

V(i,s,x,y) G [0,T] 2 x K 2 . Moreover, we have the following proposition. 
Proposition 1.2 V/i > ; 77 > 0, and Vif compact subset ofW 1 

limPfe^sup sup T^f-T— , > 77) = 0. (9) 

According to [21 Theorem 2.1.1], one can show that the limit Q holds exponentially fast as 
e -»• 0. 
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1.2 Slowly decorrelating assumption 

In this paper we are interested in a random potential with long-range correlations. Let us 
introduce some additional assumptions on the spectral gap q of the spatial power spectrum 
([3]), in order to give slowly decaying correlations to the random potential V. 

Let us note that Vi > 0, the random field V(t, •) has spatial slowly decaying correlations. 
In fact, if we freeze the temporal variable, the autocorrelation function of the random potential 
V(t, •) is given by 

i2(t,x)=E[7(t,x + y)V(t,y)] = J dp%{p)e'^ 

where Ro(p) is assumed to have a singularity in 0. Therefore, R(t, •) G" L 1 (JR ci ) and (V(t))t>o 
models a family of random fields on M. d with spatial long-range correlations which evolves with 
respect to time. However, since ([T]) is a time evolution problem, we have to take care of the 
evolution of the random perturbation V with respect to the temporal variable. In fact, as we will 
see, if the family (V(i))t>o has rapidly decaying correlations with respect to time the evolution 
problem (JTJ) behaves like in the mixing case addressed in [9J. In a time evolution problem with 
random perturbations with rapidly decaying correlations in time, even if at each fixed time the 
spatial correlations are slowly decaying, the resulting time evolution problem behaves as if it 
has mixing properties. Consequently, we have to introduce a long-range correlation assumption 
with respect to the temporal variable. Let us note that V(s,x, y) G R + x M. 2d 

[ +0 ° dt\E[V(t + s,x + y)V{s,y)]\ = +oo ^ / dp^^- = +oo. (10) 
Jo 1 J 0(pj 

Consequently, throughout this paper we say that the family (V(t))t>o of random fields with 
spatial long-range correlations has slowly decaying correlations in time if 

dp^Hl = +QO , (11) 

0(p) 

and rapidly decaying correlation in time otherwise. For technical reasons in the proof of The- 
orem 12.21 we assume that 

II I 1 2 I 1 3 

/ dpR (p)^\ + / dpR (p)^- + / dpR (p)^- < +oo. (12) 
J 0(P) J (P) J (P) 

This assumption is satisfied if |p| = C(g(p)) asp-)-0 for instance since Rq G L 1 (M d ). For the 
sake of simplicity, we assume throughout this paper that 

~ rjMi as p ^ 0, with 9 G (0, 1), (13) 

fl(p) IpI + 

and a > 0. This last assumption is used in the proof of Theorem 13. 1[ For example, in [TU] the 
authors have supposed that 

0(p)=Hp| 2/3 and R (p)= |p| ^_ 2 - (14) 

Here, a is a continuous function with compact support and a(0) > 0, v > 0, (3 G (0,1/2], 
a G (1/2, 1), and a + f3 > 1. These assumptions permit to model a random field V(t,x) with 
spatial long-range correlations for each time t > and with slowly decaying correlations in 
time. 
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1.3 Wigner transform 



To study the asymptotic phase space energy propagation of the solution <j) e of the Schrodinger 
equation d2J), let us consider the Wigner distribution of 4> e averaged with respect to the ran- 
domness of the initial data: 

W e (*,x,k) = 7 -^ 3 / dy^C)e ik - y e (t,x-e^c)^(t,x + e^c). 

(2ir) a jRd xS V 2 / v 2 / 

We refer to |23tl29| for the basic properties of the Wigner distribution. The Wigner distribution 
satisfies the following evolution equation 

W e (t,X,k)+k- V x W e (i,X,k) = 

1 r Vf^dp) . lf( p , ( pu ( 15 ) 

4^ Jri (27r) d i V e V 2/ V ' ' 2 



1 + 7 

e 2 



with initial conditions W e (0, x, k) = Wb. e (x, k), where Wo j£ is the Wigner transform of the initial 
data </>o, e of ([2]). Equation (fT5j) can be recast in the weak sense as follows: VA G Co°(M 2d ), 



/"* 1 

(W e (t), A) L2(R2d) - (W e (0), A) L2(R2d) = y (W e , k • V X A + -^£ e (s)A) L2(R2d) ds, 



where 



x, k) _ / rj (A(x, k - f) - k + f )). (16) 

Here, C^°(IR M ) stands for the set of smooth functions with compact supports. In what follows, 
we assume that Wo, e converges weakly in L 2 (M. 2d ) to a limit Wo, that is, 

VAGL 2 (IR M ), lim(Wo, £ ,A) L2(R2(i) = <W ,A) L2(R2d) . (17) 

Let us recall that by the Banach-Steinhaus Theorem Wo <e is bounded in L 2 (M. 2d ) uniformly in 
e. The main reason to introduce the additional randomness through the initial data 0o,e an d 
(S, /j,(dq)), is to make possible the weak convergence (fT7|) . Let us give an example of such a 
situation. Let /u(q) be a nonnegative rapidly decreasing function such that ||M||ii(Rd) = 1. Then, 
(M. d , f/,(dq)) is a probability space, where fJ,(dq) = /z(q)dq. Let </>o,e( x , q) = 0o(x)e -iq ' x / e , and 
0o be a sufficiently smooth function. In this case, we have || Wo^H^^d) < ||HU 2 (]R d ) II^oIIl 4 ^)) 
and Wo, e converge weakly in L 2 (M. 2d ) to Wo(x, k) = ^ydM(k)|</>o(x)| 2 . In this example, the 
initial data 0o,e(x, q) depends on the phase vector q, which is distributed according to the 
probability distribution /i(dq). Then, the averaged Wigner transform consists in taking the 
expectation of the Wigner transform with respect to the probability distribution fJ>(dq) of the 
phase vector q. 



2 Asymptotic Evolution in Long Range Random Media 

The two following sections present the asymptotic behavior of the phase and the phase space 
energy density of the solution of ([2]). The results in these two sections permit to show the 
qualitative difference between the random effects induced on a wave propagating in long-range 
random media in time and in rapidly decorrelating random media in time (|1U|) . 
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2.1 Phase Evolution 



In this section we present the asymptotic behavior of the phase of (fr t solution of ([2]). Theorem 
12. II stated below has been shown [TO] in the case 7 = 0. Nevertheless, the proof of Theorem 12, II 
remains the same as the one of |1U| Theorem 1.2], but we state this result in the case 7 € (0, 1) 
in order to provide a self-contained presentation. Under the long-range correlation assumption 
in time and medium parameters given by (|14|) . we have the following result. 

Theorem 2.1 Let us assume that the autocorrelation function R(t,x) of the random pertur- 
bation is given by ([3]), (0), and (fl%|) . and let 



a + 2/3 — 1 j kq 

' 

Let us consider the process Cc 7 ,e(*>£) defined by 

1 - / £ \ 
Cre 7 ,e(£,£) = d/(2« 7 )^ K -r> e ^' 1/(2k 7 ) J 



with 



and where <f> satisfies 



X 



: 1/(2k 7 ) ' F l/(2« 7 ) 



0(0, x) =0 Q (x). 
T/ien, /or eac/i t>0,^M d /ixed, and 7 G (0, 1), C« 7 ,e(*, 

converges in distribution to 
Co(t, = 0o(x) exp (iy/Dfa,Z)B^ (t)) , 
where (B K (t))t is a standard fractional Brownian motion with, 

a(0)Q d f+°° e~ up213 

Here, is the surface area of the unit sphere in M. d , and ei G S ' -1 . 

This result means that the phase of the field 4> evolves on the scale £ _1 /( 2k t) which is smaller 
than the one on which evolves the phase space energy density as we will see below in Theorem 
12.21 The random perturbations induce a random phase modulation of the fields <p given by a 
fractional Brownian motion. This result shows an important difference with the case of media 
with rapidly decorrelating random perturbations in time (|10p . for which the phase and the 
phase space energy density evolve on the same scale e" 1 [TO]. 

Let us note that the scale of evolution of the random phase is shorter in the case 7 > 
than in the case 7 = 0, that is e -1 /^ 7 - 1 <C e~ 1 /( 2re °) -C e _1 . This fact comes from the random 
perturbations which are stronger in the case 7 > than in the case 7 = 0. These strong 
perturbations in time lead to the development of a random phase modulation earlier than in 
the case 7 = 0. However, the scale of evolution of the phase space energy distribution is the 
same e _1 in both cases, 7 > and 7 = 0. 
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2.2 Limit Theorem 



This section presents the main result of this paper. Theorem 12.21 describes the asymptotic be- 
havior of the phase space energy distribution of 4> t solution of the random Schrodinger equation 
d2]). First, let us introduce some notations. Thanks to (fT7|) . let r = sup e || Wo,e|lL 2 (R 2d ) < °°> 



B r = {A G L 2 (R 2d ),\\X\\ L 2 



< r 



be the closed ball with radius r, and {g n ,n > 1} be a dense subset of B r . We equip B r with 
the distance djg r defined by 

+oo j 
3=1 

V(A, /i) G (£v) 2 > so that (B r ,dQ r ) is a compact metric space. Therefore, (W e ) e is a family of 
process with values in (B r ,djs r ), since ||W / e(0llL 2 (R 2d ) = ll^o.elliZf 



Theorem 2.2 T/ie family (W e ) (Eg (o,i); solution of the transport equation (|15p . converges in 
probability on C([0, +oo), (£> r , ds r )) as e — >■ to a Zz'mz'i denoted by W. More precisely, VT > 
and \/rj > 0, 



limP sup d Br {W e {t),W(t)) > 7] = 0. 
\te[o,T] / 

W is the unique weak solution uniformly bounded in L 2 (M. 2d ) of the radiative transfer equation 

d t W + k • X7 X W = CW, (18) 
with W(0,x, k) = Wo(x, k). Here, C is defined by 

CcpQs.) = J dpa(p - k) ( v (p) - ^(k)) , (19) 

with <p £ C°°(R d ) and where 

a (p) = 2 ^o(P) . (20) 

(27r) d B (p) V j 

Here, Co°(lR rf ) stands for the infinite differentiable functions with compact support. Let us 
precise that W is a weak solution of the transfer equation (|18|) means that 



(W(t), X) L2(Rd) - {W , X) L2(Rd) - J n (W(s),k • V X A + CX) L2(Rd) = 0, 
Vt > and VA G Cg°( 



The topology gererated by the metric dg r is equivalent to the weak topology on L 2 (M. 2d ) 
restricted to B r . We cannot expect a convergence on L 2 (R. 2d ) equipped with the strong topology. 
In fact, the following conservation relation [|W e (f)||i2mad> = [|WQ )6 [|i2ma<i> cannot be satisfied 
by the limit W. Moreover, let us note that the Wigner distribution W £ is self-averaging as e 
goes to 0, that is the limit W is not random anymore. This self-averaging phenomenon of the 
Wigner distribution has already been observed in several studies [B] and is very useful in 
applications. 

The proof of Theorem [221 is given in Section [B] and is based on an asymptotic analysis using 
perturbed-test-function and martingale techniques. 

Equation (|18|) describes the asymptotic evolution of the phase space energy distribution of 
the field cp e solution of the random Schrodinger equation ^fy. The nonlocal transfer operator 
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C defined by (fT§|) describes the energy diffusion caused by the random perturbations, and the 
transfer coefficient o~(p — k) describes the energy transfer between the modes k and p. 

Let us note that the result of Theorem 12.21 does not depend on whether J dpc(p) is finite 
or not. In other words, the radiative transfer equation (|18|) is valid in the two case, slowly 
and rapidly decaying correlations in time (fTUj) . Consequently, the phase space density energy 
distribution evolves on the same scale e -1 in both case, rapidly and slowly decaying correlations 
in time. However, as we will see in Section [31 the solutions of these equations in the two cases 
behave in different ways. As it has been discussed in Section fl. 21 in the case of rapidly decaying 
correlations in time (|iup . that is J dpcr(p) < +oo, the radiative transfer equation (|18|) has the 
same properties as in the mixing case addressed in [S] . In the case of slowly decaying correlations 
in time (|l(jp . that is J dpcr(p) = +oo, we observe a regularizing effect of the solution of (|18|) 
which cannot be observed in the case of rapidly decaying correlations in time. As we will see 
in Theorem 13. 11 the unique solution W of (fTHj) is actually the unique classical solution of ([T5]). 



3 Regularizing Effects of the Radiative Transfer Equation ( TTSj ) 



In this section we investigate the regularizing effect of the radiative transfer equation (|18|) . We 
show that the solution of (|18|) becomes instantaneously a smooth function, that is upon t > 0, 
despite the nonsmoothness of its initial data. 

Theorem 3.1 Let W G L 2 (R 2d ) and W be the unique weak solution of the radiative transfer 
equation (|18|) . Then, under f|13j) and Vio > 0, we have 



WGC°((0,+oo), f| H k (R 2d )) nL°°([t ,+oo), f| H k 



k>0 k>0 

so that we also have 

W G C°((0,+oo),C°°(IR M )) nL oo ([t ,+oo),C oo (M M )). 
Moreover, W is given by 



where 



*(q) = J dp<7(p)(e*-«-l). (21) 



In TheoremlXTl H k (R 2d ) stands for the kth Sobolev sapceonR . The proof of Theorem I3TT1 
is given in Section [3 This result is an important consequence of the slowly decaying correlation 
assumption in time (|11|) and cannot be observed in the case of rapidly decaying correlations in 
time (|l(jp . Theorem 12.21 and Theorem 13. H leads us to the following corollary. 

Corollary 3.1 The family (W e ) eg (o,l)j solution of the transport equation (|15j) with Wq G 
L 2 (R 2d ), converges in probability on C([0, +oo), (B r , dfl r )) as e — > to a limit denoted by W, 
which is the unique classical solution uniformly bounded in L 2 (M. 2d ) of the radiative transfer 
equation 

d t W + k • V X W = CW, (22) 
with W(0, x, k) = Wo( x ik), and where C is defined by (fT9|) . 

Thanks to the long-range decorrelation assumption in time (|13|) . the limit obtained in 
Theorem 12.21 in not only a weak solution of the radiative transfer equation (|18p but also a 
classical solution. In Section HI we give an explanation of this regularizing effect using a 
probabilistic representation of the solution of (|18|) in term of a Levy process with jump measure 
a(p)dp. 
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4 Probabilistic Representation 



In this section we discuss the probabilistic representation of the solution of the radiative trans- 
fer equation (|18|) in term of Levy process. A Levy process is a stochastic process with inde- 
pendent and stationary increments. We refer to [U [33] for the basic properties of the Levy 
processes. A particular property of these processes is that they are entirely characterized by 
their characteristic exponent defined by E[e jq ' Lt ] = e - **( q ). For instance, there exists a Levy 
process associated to the generator (fH?|) and for which its characteristic exponent is given by 



Proposition 4.1 Let W be the unique weak solution of the radiative transport equation (|18p 
with initial datum Wo G L 2 (M. 2d ). Then, there exists a Levy process (Lt)t>o with caracteristic 
exponent (|2ip such that Lq = 0, and 

W(t, x, k) = E [W Q (x - tk - jf* L s ds, k + Lt) . 

According to Theorem \S.l\. W is the unique classical solution uniformly bounded in L 2 (M. 2d ) of 



The proof of Proposition 14.11 is given in Section [5] According to ([2"T|) . the Levy process 
(Lt)t is a pure jump process with jump measure <r(p)dp. Therefore, because of the long-range 
correlation property (|13p the jump process {Lt)t has infinitely many small jumps. This last 
property of the symmetric process (Lt)t is the key to show the regularizing effect of the radiative 
transfer equation (fT8"|) . In fact, according to [Ml Proposition 1.1] L t has a smooth bounded 
density, which permits to obtain smoothness to the variable k. Moreover, the transport term 
in (|18|) permits to transfer the smoothness of the variable k to the variable x. The probabilistic 
representation of the solution of ([15]) presented in Proposition 14.11 will be usefull in the proof 
of Theorem 15.11 



5 Fractional Radiative Transfer Equation 



The goal of this section is to give an approximation of the radiative transfer equation (|18p . 
This approximation permits to get a simpler radiative transfer model and derive an explicit 
formula of the solution, from which one can easily extract the exponent of the decaying power 
law. 

The nonlocal effect of £, defined by (|19|) . is characterized by the support of the transfer 
coefficient c(p) defined by (|2U|) . The extension of this support is tantamount to the extension 
of the nonlocal effects of the transfer operator C. In this section, we are interested in what does 
the radiative transfer equation (|18|) look like if the support of the transfer coefficient increase? 

Let us recall that the transfer coefficient is given by 



<np) 



fl(p) 



12(0, p) 



E[V(t + s,x + y)V(s, y)] e - ip - x dxdt, 



which is the power spectrum of the two-point correlation of the random potential V at frequency 
0. In order to study the impact of the extension of the support of <r(p), let us consider the 
power spectrum given by 



1 



E 



t + s x + y^ yf* y 



rj r\ 



e- ip - x dxdt = V d+e a{i]p). 



(23) 



Here, the parameter rj represents the extension of the support of a v (p), and therefore the 
extension of the nonlocal effects of the corresponding transfer operator jC 11 as rj goes to 0. Let 
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us remark that this scaling corresponds to the long space and time diffusion approximation of 
the radiative transfer equation (fTgj) . 

In this paper, we have assumed (|13p . which means that that the transfer coefficient behaves 
locally at like the transfer coefficient of the fractional Laplacian — (— A) e / 2 . Thanks to the 
scaling (|23p this local behavior becomes global as shown in Theorem 15.11 In fact, using (|13p 
we have 

Vpel^\{0}, lim^( P ) = ^^. (24) 
Moreover, let us assuming that Rq is bounded on each compact of R d \ {0}, we also have 

< <7"(p) < VpGM d \{0}, (25) 

with C > 0. 

Theorem 5.1 Let us assume that Rq is bounded on each compact ofM. d \ {0}. Then, we have 

W G L 2 {R 2d ) and \ft > 0, lim W v (t, x, k) = W°°(t, x, k) 

rj— >0 

pointwise and weakly in L 2 (R 2d ). Here, W v is the unique classical solution uniformly bounded in 
L 2 (M. 2d ) of the radiative transfer equation (|18p . and W°° is the unique weak solution uniformly 
bounded in L 2 (M. 2d ) of the following radiative transfer equation 

dtW^ + k- V X W°° = -a(e)(-A k ) e / 2 W°°, VF°°(0, x, k) = W (x, k). (26) 

Here, (— A^) 9 / 2 is the fractional Laplacian with Hurst index G (0, 1), and 

a6T(l - 6) f g 

with ei G S d_1 and T(z) = J +o ° t^e^dt. Moreover, W°° is given by the following formula 
W°°(t,x,k) = j^- d J dyd^y+^h-^I^^^Woiy^ + ty). (27) 

The proof of Theorem 15.11 is given in Section [U] and uses the probabilistic representation 
obtained in Proposition 14.11 

Equation (|26p is an approximation of (|18|) when we extend the nonlocal effect, and it 
permits to get an explicit formula (|27|) of its solutions. This formula permits to exhibit a 
damping term, which means that the phase space energy distribution decays with respect to 
time. Moreover this damping term obeying to a power law with exponent 9 G (0, 1). In wave 
propagation such a kind of result result has been already obtained for the wave equation in 
one dimensional propagation medium [22], and physical models involving fractional Laplacian 
have been proposed to predict such a power law |15j . 

Regarding the regularity of W°° we can prove the following proposition in the same way 
as Theorem 13.11 This result states that the unique solution of ([26]) W°° is actually a smooth 
function which is the unique classical solution of (|26|) . 

Proposition 5.1 Let W G L 2 (R 2d ) and W°° be the unique solution of the radiative transfer 
equation ([25]) . Then, Vto > 0, we have 

W°° GC°((0,+cx)), p| H k (R 2d )) nL°°([t ,+cx)), p| H k (R 2d )), 

k>0 k>0 
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so that we also have 

W°° G C°((0,+oo),C oo (R 2d )) nL~([t ,+oo),C°°(K 2d )). 

Then, W°° is the unique classical solution uniformly bounded in L 2 (M. 2d ) of the radiative trans- 
fer equation 

dtW^ + k- V X W°° = -a(9){-A) e / 2 W°°, lU°°(0,x,k) = W (x,k). 
6 Proof of Theorem 12.2 

The proof of Theorem 12.21 is based on the perturbed-test-function approach. Using the notion 
of a pseudogenerator, we prove tightness and characterize all subsequence limits. 

Using a particular tightness criteria, we prove the tightness of the family (W / e) e e(o,i) on the 
polish space C([0, +oo), (£> r ,dg r )). In the next section, we characterize all subsequence limits 
as weak solutions of a well-posed radiative transfer equation. 

We have the following version of the Arzela-Ascoli theorem 125] for processes with values 
in a complete separable metric space. 

Theorem 6.1 A set B C C([0, +oo), (B r ,dj3 r )) has a compact closure if and only if 

VT > 0, lim sup rriTig, fj) = 0, 

with 

m T (g,v) = SU P d Br {g{s),g{t)). 
(M)e[o,T] 2 

\t-8\<1) 

From this result, we obtain the classical tightness criterion. 

Theorem 6.2 A family of probability measure (P e ) eg / 01 ^ on C([0, +oo), (B r , dj3 r )) is tight if 
and only if 

VT > 0,r/ > lim sup F e (g; m T {g,v) > H) = °- 

'^ eG(0,l) 

From the definition of the metric dg r , the tightness criterion becomes the following. 

Theorem 6.3 A family of processes (^ £ ) e g(o,i) is tight on C([0, +00), (B r , ds r )) if and only if 
the process ({X e , X) L2(RdxRd) ) ee(0A) is tight on C([0,+oo),M) VA G L 2 (R 2d ). 

This last theorem looks like the tightness criterion of Mitoma and Fouque [63\ I2T] . 

For any A G L 2 (R 2d ), we set W £) a(*) = (W t (t), X) L2(R2dy According to Theorem E3 
the family (W t ) e is tight on C([0, +00), (£> r ,cfe r )) if and only if the family (W et \) e is tight on 
C([0, +00), R), VA G L 2 (M. 2d ). Furthermore, (W e ) € is a family of continuous processes. Then, 
according to [11, Theorem 13.4], it is sufficient to prove that, VA G L 2 (R 2d ), (W 6: \) e is tight 
on T>([0, +00), K), which is the set of cad-lag functions with values in R. Finally, using that 
the process W e is a process with values in B r and that the set of all smooth functions with 
compact support in R 2d is dense in L 2 (K. 2d ), it is sufficient to show that (W £: \) e is tight on 
P([0,+oo),R) VA G C^(R 2d ). 
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6.1 Pseudogenerator 

We recall the techniques developed by Kurtz and Kushner [28]. Let A4 e be the set of all T^- 
measurable functions f(t) for which sup t < r E [|/(£)|] < +00 and where T > is fixed. Here, 
J-f = J~ t / e and {J-t) is defined by ©. The p — lim and the pseudogenerator are defined as 
follows. Let / and f s in M e V<5 > 0. We say that f = p - lim^ f s if 

supEfl/" 5 ^)!] < +00 and limE[|/ 5 (i) - f(t)\] = Vt. 
t,s <*->o 

The domain of A e is denoted by V (A e ). We say that f £ V (A e ) and A e f = g if / and g are 
in V (A £ ) and 

\m [f(t+s)]-f(t) g{t) 



■p — lim 

<5^0 



0, 



where E^ is the conditional expectation given T\. A useful result about pseudogenerator A e is 
given by the following theorem. 

Theorem 6.4 Let f eV (A e ). Then 



is an (J 7 /) -martingale. 



M}(t) = f(t) - /(0) - J A e f(u)du 



6.2 Tightness 

Proposition 6.1 VA G C^(R 2d ), the family (W e ,A)ee(o,i) is tight on V ([0, +00), R). 

Proof According to [28, Theorem 4] and because W e is a process with values in B r , we need to 
show Lemma [6. 21 and Lemma [6.41 b elow . Throughout the proof Propostion l6.14 let A G Co°(M? d ), 
f be a bounded smooth function, and = f(W e< \(t)). 



Lemma 6.1 VT > 0, 



e[ sup WCeH^Whm^d)] < +00. 

te[o,T] v ; 



Proof (of Lemma 16. 1J) First, 

/■V2 

>.a®2d\ < /„ ,j 

'-1/2 

Let us fixe x, k, and u. Let 

<£i,A,x,k, u (p) = • V k A(x, k + up). 



£ e A(t)||£ a(RM) < I dxdkj I y(-A- )( ip) e J P x / £ p. V k A(x,k + np) 



(28) 



It is clear that 4>ix^,Ku ^ Consequently, Vj = (V, 4>\^^ u ) w H is centered Gaussian process 
with a pseudo-metric m\ on [0, T] given by 



mi(t, s) = E 



:l+7 



1/2 



Then, 



m\(t,s) < 2sup|V k A(x,k) 

x,k 



^(PIOIPIIPI 2 ) 1 ^ V(t, S )G[0,T] 2 , 
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diam 2 m ([0,T])<2 / # (p)|p| 2 | V k A(x, k + np)| 2 . 



and 



Here, diam mi ([0,T]) stands for the diameter of [0,T] under the pseudo-metric mi. According 
to [21 Theorem 2.1.3], we have 



E 







2 


sup 






t£[0,T] 





( rdiam mi ([0,T])/2 > 

<KU H 1/2 (r)dr 



<Ci 



6»i(x,k,n) / / rp 



where H(r) = ln(iV(r)), and N(r) denotes the smallest number of balls, for the pseudo-metric 
mi, with radius r to cover [0,T]. Moreover, 8\ is given by 



# 2 (x,k,n) = 2 / i? (p)|p| 2 |V k A(x,k + np)| 2 . 



Consequently, 



i(x,k,n) / / rp 

dx.dk | / Win I C 2 2 1+ ] 1 < +oo, 



o 



since Ro and A have a support included in a compact set, and that concludes the proof of this 
lemma. □ 

Thanks to Lemma 16.11 we can use the notion of pseudogenerator introduced in Section 16.11 
Therefore, we have 



= nw € , x (t)) 

where C t is defined by ([To]) , and 



e 2 



L 2 



Ai(x,k) = k- V x A(x,k). 



(29) 



Let 



1 



1 + 7 

e 2 



/'(W e , A (t)) / dxdkW £ (t,x,k) 



" + 00 



^(^p) 



(2vr) d i 

Lemma 6.2 VT > 0, and n > 



,i(-u-t)p-k/e ip-x/e 



(A(x,k 



f)-A(x, k + |; 



limP( sup |/! e (t)| > r?) = 0, and supE [|/i e (t)|] = 0(e (1 " 7)/2 ). 

e V 0<i<T 7 t>0 



(in. 



Proof (of Lemma 16. 2|) Using we have 



1 + 7 



with 



fl(t) = —f\w^{t))(w e (t),c u \{t)) L2{M:2iy 



& e A(t, x, k) = ^i— / _^Mn^PL e *P-x/e A ( k _ £) _ A( k + £ 
' V ; (2vr) d U g(p) - ie-yk ■ p \ 2 7 v ' 2' 
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Lemma 6.3 

E 



sup ||£i i£ A(t)|| L 2 
te[o,T] 



1 + 7 



sup ||£i i£ A(t)||| 2 
te[o,T] 



< +00 and lime 2 E 

e 

Proof (of Lemma I6.3j) Here, we use the same argument as in Lemma 16.11 First, 



0. 



Let us fixe x, k, and u. Let 

02,A,x,k,w(p) 



Jp-x/e 



j(p) - ie^k • p 



p • VkA(x, k + up). 



(30) 



According to (fT2"j) </>2.x,k,u £ Consequently, V2 = (V, 4>2,x,x,k,u) w ^ is centered Gaussian 
process with a pseudo- metric 777-2 on [0, T] given by 



m2 (i, s) = E 



(wfe^-wb*))' 



1/2 



Then, 



m|(t,«) < 2sup|V k A(x,k)| 2 / d P R (p 



x,k 



,0(p)Ip| 2 \ 



2 (P) ) ^ 



V(i,s) e [0,7f 



and 



Pi 



&m; 2 ([0,T]) <2 / i?o(p)^T|V k A(x,k + up)P 

2 (p) 



Here, e?iam m2 ([0, T]) stands for the diameter of [0, T] under the pseudo-metric 7772. According 
to [21 Theorem 2.1.3], we have 



E 



sup 

te[o,T] 



<Ci 



where 



Consequently, 



e 2 (x,k,-u) 



P 



'In C- 



T 

2 r 2 e l+7 



^(x,k,u) =2 / J R (p)^-|V k A(x,k + 7 J p)|^. 



2 (PJ 



dx dp 



e 2 (x,k,«) 







'In C 2 - 



dr < +00, 



r 2 e l+7 

since Rq and A have a support included in a compact set, and that concludes the proof of 
Lemma 16. 31 D 

Then, the proof of Lemma 16.21 is a direct consequence of Lemma I6.31 D 

The following lemma is a consequence of Lemma 16.31 Lemma 16.51 and Lemma [ST 
Lemma 6.4 VT > 0, {A e (/q + ft) (*), e 6 (0, 1), < f < T} is uniformly integrable. 



Proof (of Lemma 16. 4|) First, let us show that we can compute the pseudogenerator at /0+/1 • 
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Lemma 6.5 

sup E\\\£ e (£ u \(t))(t)\\ 2 L2 fR2d) < + x. 

65(0,1) L 

Proof (of Lemma 16.5 j) We have 
£ e (£i, e A(i))(i,x,k) = 

y(^,d Pl )y(- T ^, ( ip 2 ) e ^+^)- x A 



.(A(x,k-£-?)-A(x,k-£ + ?)) 



S (p 2 )-^(k-^)- P2 v v ' 2 2 ; v ' 2 2 

(A(x,k + ^-^)-A(x,k + ^ + ^)) 



(p 2 )-i e 7(k+^) .p 2 

Let us note that 

E[V(h, d Pl )V(t 2 , d P2 )V*(t 3 , dp 3 )V*(U, dp 4 )] = 

(27r) 2d i?(ii - t 2 , pi)£(i 3 - U, p 3 )*(pi + P 2 )5(p 3 + p 4 ) 
+ (2n) 2d R(h - t 3 ,pi)-R(t 2 - i4,Ps)5(pi - p 3 )5(p 2 - p 4 ) 
+ (2TT) 2d R(h - t 4 , pi)A(t 2 - i 3 , Ps)*(pi - p 4 )<5( P2 - P3 ). 

Then, using the smoothness of A and (|12|) . we obtain 

E[[|£ e (£i, e A(i))(i)|| 2 2(R2d) ] <c|y d Pl dp 2J Ro(pi)i?o(p 2 ) 

( l Pl L sup|V k A(x, k)| + ^ Pl ^ sup ||-D 2 A(x, k 
A 0(Pi) x,k fllPlJ x,k 



(31) 



1 2 

x sup| V k A(x, k)| + -i^L sup || jD 2 A(X) k) 

V 0(P2j^ x,k 0(P2j x,k 

+ ( IP 2 ' P2 f sup|V x A(x,k)| + -J— (| Pl || P2 | + |p 2 | 2 )sup||/J 2 A(x,k)||V 

V (P2j x,k 0(P2) x,k ' 

+ ( IP 2 ' P2 f sup|V x A(x,k)| + -J_(| Pl || P2 | + |p 2 | 2 )sup||D 2 A(x,k) 

V (P2j x,k 0(P2j x,k 

x ( |P2llPi| 2 sup | VxA(X)k )| + _J_(| p2 || pi | + | Pl | 2 ) SU p||D 2 A(x,k 

V (Pi) X,k fl(Pl) x,k 



< +00, 



since A and i?o have a support included in a compact set. That concludes the proof of Lemma 
□. 



Lemma 6.6 

sup E[||£ e A(t)|| 2 2(R2d) x ||£ lie A(t)|| 2 2(R2d) ] < +00. 



Proof (of Lemma 16. 6|) This result follows from the temporal stationarity of the process 
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(V(t,-)) t . We have 

E[||£ e A(t)||^ 2 ( R d xR d) x ||£i, £ A(t)||^ 2 ( Rd><Rd) 
< J J dxidki(ix2C?k2 J 



1/2 ,1/2 

du\ I du2 

/2 7-1/2 

E[|(y(Ve 1+7 ),0 1 ,A ) x lj k ll u 1 > w ,^| 2 Ki>(Ve 1+7 )>2Ax 2l k 2! « 2 >^| 2 ] 

/■/■ ri/2 fi/2 

< / / dxi<iki(ix2dk2 / dui / du2 

JJ J -1/2 7-1/2 

where <^i )Al x,k,u and <^2,A,x,k,u are defined respectively by ([25]) and ([3D]) . Moreover, 
E[|(^(^ 1+7 )^l,A,x ll k 1 ,« 1 >^^| 4 ] =3E[|(y(0)^ l!A)X , k)M > w ^| 2 ] 2 



< 3 



^o(p)|p| 2 |V k A(x,k + «p) 



< +oo, 



and thanks to (fT2j) 



E[|<y(Ve 1+7 ),0 2 ,A,x 1 ,k llUl ) w , iW | 4 ] =3E[|(t/(O),0 2 , A , x , M ) w/jW | 2 ] 2 



< 3 



^o(p) 



2 (p) 



|V k A(x, k + up) 



i 2 



< +oo. 



That conclude the proof of Lemma 16.61 since A and i?o have a support included in a compact 
set.D 

Consequently, thanks to Lemma 16.3} Lemma 16.5} and Lemma 16.61 we have 
■Wo + /!)(*) =f'(W e , X (t)) [W eM {t) + {WS),Ce{Cl,e\{t)){t)) L2{v?d) 

+ a 2 /(Ty £ , A (t))(^(t),A,eA(t)) L2(R2d) <^(t),£ e A(t)) i2(R2d) 

+ 0( e (l+7)/2), 

and sup e f E[|^4(/g +/f)(i)| 2 ] < +oo. That conclude the proof of Lemma and then the proof 
of Proposition 16.11 □ ■ 

6.3 Identification of all subsequence limits 

To identify all the subsequence limits of the process (W e ) e , we show that all such limit processes 
are a weak solution of a deterministic radiative transfer equation. Let us note that in this case 
all the limit processes are deterministic. This fact means that the convergence also holds in 
probability. In the next section, we will see that this transport equation is well posed. In 
particular, this will imply the convergence of the process (W e ) e himself to the unique solution 
of the radiative transfer equation. 

Proposition 6.2 Let W be a accumulation point of (W e ) e . Then, W is a weak solution of the 
radiative transfer equation 

d t W + k • V X I^ = CW, 
with W(0, x, k) = Wo( x )k). Here, C is defined by 



Op(k) = J dpa(p - k)(y>(p) - <p(k)) = J dp<7(p)(p(k + p) - y>(k)) Vy> G C ° 
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with 



2R (p) 



(27r) rf 0(p)' 

Proof (of Proposition 16. 2|) In this proof we use the following notation 

(p ® V( x i> k i> x 2,k 2 ) = v?(xi,ki)^(x 2 ,k 2 ). 



Let 



where 

ffi, e (t,x,k) 
1 



(27r) M i 2 



=^/(^, A (t))<We(t), ffl, e (t)) i2(R2d) 

+ a$f(w etX (t))(w e (t) ® w e (t), # 2 , e (t)> L2(R4d) , 



'du // (E t e [t>(u/e 1+7 ,dpi)1>(Ve 1+7 ^P2)] -E[y(0,dpi)y(0,dp 2 ) 



x e «(pi+P2)-x/e e i(u-t)(p 1 +p 2 )-k/e 



1 



(p 2 )-ie7(k-^)-p 2 
1 



v ' 2 2 1 v 2 2' 

Pi P2 



g(p 2 ) -ie^(k + 



: f A (x,k+ El _ ») - A(x,k + HI + E 2 . 

k + ^-p 2 ^ V 2 2 ; v ' 2 2 



and 



i? 2j£ (t,xi,ki,x 2 ,k 2 

"+00 



(2vr 



Wi ( E t[^Ke 1+7 ,dPl)^(«A 1+ ^dp 2 )]-E[t/(0,d Pl )F(0,dp 2 )]) 



Q ipi-xi/e „ip 2 -x 2 /e i)(pi-ki+p 2 -k 2 )/e_ 



1 



fl(pi) - ^e 7 ki • pi 



x(A(x 1 ,k 1 -^)-A(x 1 ,k 1 + ^; 
x(A(x 2 ,k 2 -^)-A(x 2 ,k 2 + ^)). 



However, according to ^ 
t 



E 



and 



^(«+3F. d Pi)V"(«+3F» d Pa 



E 



-(fl(px)+fl(pa)) 



y(o,d Pl )y(o,dp 2 ) 
t 

7 



U ^(^,dpi)^(^,<W) - (2vr) rf e~ 2 0(P 1 )«l?o(pi)5(pi + P2 ), 



E 



y(o,d Pl )y*(o,d P2 ) 



e -( fl (Pi)+fl(P2))uy(_ )dpi )y*(__ )dp2 ) _ (2vr)V 2 ^)«l?o(pi)5(pi - P2 ) 



Consequently, 



/!(*) =e 1+7 [^/(Ty ejA (t))(W e (t),F l!e (t)> i2(R2(i) 

+ a*/(W e , A (t))<W' e (t) ® W e (t), H 2)e (t)> i2(K4rf) 
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where 

#L, 6 (t,X,k) 



(27r) M i 2 



V(t/e 1 ^,dp 1 )V(t/e^,dp 2 ) 



1+7 



i(pi+p 2 )-x/e 



(fl(Pi) + fl(P2) - i^k ■ (pi + p 2 )) 



1 



(p 2 )-ie7(k-PA).p 2 
1 



A(x,k 



Pi P2 



0(p 2 )-ieT(k+^) -p 2 



+ 



(2vr) c 



c/p 



#o(p) 



2g(p) L fl (p) + t^(k-f).p 

/ x , ■ -v/, i ( A ( X > k + P) - A ( X 

fl(p) +«eT(k + f) p V 



x,k) - A(x,k - p) 



and 



#2 ! e(* ) X 1 ,ki,X2,k 2 ) 



(2vr) M i 2 



V^t/e^dpiMt/e^cipa 

e ipi-xi/e e ip2-X2/e 



(g(pi) - ie^ki • pi)(g(pi) + fl(p a ) - ie 7 (ki • pi + k 2 • p 2 )) 



x(A(x 1 ,k 1 -^)-A(x 1 ,k 1 + ^; 
x(A(x 2 ,k 2 -^)-A(x 2 ,k 2 + ^; 



1 



+ / d P 



^ (p)e ip - (xi_X2)/e 



(2vr) d 7 ^ (g(p) - *£7ki • p)(2 (p) - *e7(ki - k 2 ) • p) 



A(xi, ki - -) - A(xi, ki + -) 
A(x 2 , k 2 + -) - A(x 2 , k 2 - -) 



Lemma 6.7 



supE[|/J(t)|]=0(e 

t>o 



l+7> 



Proof (of Lemma 16. 7|) In the same way as in Lemma l6.5( using (|12]) . ()3ip . the smoothness 
of A, and that A and Rq have a support included in a compact set, we have 



supE[||#i £ 



< C / / dpidp 2 i?o(pi)^o(p 2 ) 



( l Pl L sup| V k A(x, k)| + sup ||£> 2 A(x, k) 

v 0(PiT x,k (Pi) x,k 

x ( l P /^' \ sup|V k A(x, k)| + ^ r sup ||£> 2 A(x, k 



+ 



IP1IIP2I 



0(Pl)0 2 (P 2 ) x,k 



sup|V x A(x,k)| + 



2 (P2j x,k 

PlHP 2 | 



+ 



|P2| 



sup||D 2 A(x,k)|| 



+ ( Jp i 1 |feL\ ^P|V x A(x,k)| + 



L)0 2 (P2) x,k 
|P2||Pl| 2 

0(p2)g 2 (pi) xTk 



fl(Pl)fl(P2) 2 (P2)^ x,k 

2 

sup||L> 2 A(x,k)|| 



IP1IIP2I |P2| 



fl(Pl)fl(P2) 2 (P2)^ x,k 



sup|V x A(x,k)| + 



IP2IIP1I iPil 



fl(P2)fl(Pl) 2 (Pl)^ x,k 



sup||L> 2 A(x,k) 



< +00, 



and 



supE[||tf 2 , e 



< C // dpidp 2J R (pi) J Ro(P2] 



|Pl| 2 |P2| 2 



2 (Pl)0 2 (P2) x,k 



sup|V k A(x,k)| 4 . 
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That concludes the proof of Lemma 16.71 □ 

Consequently, after a long but straightforward computation, we get 

A\r Q + ft + fl)(t) =d v f(W e!X (t))[W eM (t) + (W e (t),G lte X) L2{R2d) } 

+ d 2 J(W e , x (t))(W e (t) ® W e (t), G 2 , e A) L2(R4d) (32) 

+ 0(e (l+7)/2 ); 

where 

Gi,A(x, k)=-^] d P R ( P ) [ 0(p) _ Z£7(k _ f) . p (A(x, k) - A(x, k - P) ; 

' / A(x,k + p)-A(x,k) 



fl(p)-*^(k + |)-p 
1 /• g(p)^ (p)e ip(xi " X2)/e 



G a , e A(x 1 ,k 1 ,x 3l k B ) =- j^J rfP (fl(p) _^ kl . p)(2fl(p) _ ie7(kl _ k2) . p) 



x (A(x 1 ,k 1 -|)-A(x 1 ,k 1 + |; 
x (A(x 2 ,k 2 + |) -A(x 2 ,k 2 -|; 



and 



G 3>e A(x 1 ,k 1 ,x 2) k 2 )=-^ I j d P( (p)_^ kl .p)(2 (p)-, e 7( kl - k2 ). p 

x (A(x 1 ,k 1 -|)-A(x 1 ,k 1 + |) 



x 



A(x 2) k 2 + |)-A(x 2) k 2 -|)). 



Lemma 6.8 

lim ||G 2ie A|| L 2( R 4d) + ll^eAll^a^d) = 0. 

Proof (of Lemma I6.8j) We will just show the result for G 2j£ , the proof is exactly the same 
forG 3 , e . 

||G 2 , e A||2 2(R4d) = J dx 1( ix 2 d kl dk 2 | J dp^(x 1 ,x 2 ,k 1 ,k 2 ,p)e i P ( Xl - X2 )/ £ 

where 

g(p)-Ro(p) 



* e (xi,X2,ki,k 2 ,p) 



(g(p) - ie^k! • p)(2g(p) - »£7(ki - k 2 ) • p) 

rl/2 rl/2 

xp - / duiVkA(xi,ki + uip) x p • / du 2 V k A(x 2 , k 2 + u 2 p). 

7-1/2 J-l/2 



'-1/2 J-l/2 

Let us recall that i?o an d A have their supports included in a compact set. Then, 

Mp)\p\ 2 



with 



sup |* e (xi,x 2 ,ki,k 2 ,p) - ^ (xi,x 2 ,ki,k 2 ,p)| < e 7 G- 
xi,ki,x 2) k 2 01PJ 



#o(p) f 1/2 f 1/2 
* (xi,x 2 ,k 1 ,k 2 ,p) = — — — p • / ditiV k A(xi,ki + uip) x p • / du 2 V k A(x 2 , k 2 + u 2 p) 

2g p J -in J -in 



2g(p) J J -i/2 
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Moreover, for almost every (xi, x 2 , ki, k 2 ) 

lim J dp*(x 1 ,x 2 ,k 1 ,k 2 ,p)e ip ( Xl - X2 )/ £ = 

by the Riemann-Lebesgue lemma, since i?(p)^y G Consequently, the result of Lemma 

16.81 follows from the dominated convergence theorem. □ 



Lemma 6.9 We have 
with 



lim\\(G l , e -G 1 )X\\i m2d) =0, 



G x A(x, k) = -i- / dp^- (A(x, k + p) + A(x, k - p) - 2A(x, k) 



' /dp^^(A(x,k + p)-A(x,k) 
J P v 



(2ir) d J * ( p ) 

^ I 1 2 

Proof (of Lemma 16. 9|) This Lemma is a direct consequence of R(p) f) \ G L 1 (M d ). □ 



Let = /o(i) + /i (*) + /!(*)■ By Theorem^! (M| E (i)) t>0 is an (J^-martingale. That 
is, for every bounded continuous function $, every sequence < s% < ■ ■ ■ < s n < s < t, and 
every family {fJ-j)j£{i,... >n } £ L 2 (U. d ) n , we have 



E 



HW e>lij ( Sj ),l < j <n)(r(t) - f(s) - f A e f(u)du) 



0. 



Using (|32p . Lemma l6.2| Lemma l6.7| Lemma l6.8| and Lemma 16.91 we obtain that 
M f)X (t) =f(W x {t)) - f(W x (0)) 

- /\/(W A («))[W Al («) + (^(n),G 1 A) i2(R2d) ] 

«* 

is a martingale, and where Ai is defined by (|29|) . More particularly, Let us consider / be a 
smooth function so that f(v) = v, Vu such that \v\ < r\\X\\ L 2m2d\. Then, 



M\(t) =W x (t) - W x (0) 

du[W Xl (u) + <W»,GiA} L2(R2(i) ] 

is a martingale with a quadratic variation equal to 0. Consequently, Mx = 0, that concludes 
the proof a Proposition 16.21 ■ 



6.4 Proof of the weak uniqueness of the radiative transfer equation (fTE)) 

In this section we show the uniqueness of weak solutions of the radiative transfer equation (|18p . 
To this end will construct a good test function using an approximation Cm defined 

C N (p{k) = dpo-(p)(<p(p + k)- <p(k)). 

J\p\>l/N 

of C defined by 

Ctp(k) = J dpa(p) (<p(p + k) - <p(k)) , 

V99 G Cq°(K ). £ at is an approximation of C with a finite scattering coefficient Jjpi>i/jv dpa{p). 
As we will see Cm is a well suited approximation of C to construct a test function which enables 
us to show the weak uniqueness of the transfer equation. 
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Proposition 6.3 Let A G Cg°(R M ). ViV > 1, i/iere exists a unique solution Wjv of i/ie radiative 
transfer equation 

d t W N (t,xi,k) =k- V x W JV (t ! x,k)+£ iv M A Jv(t,x ! k) (33) 

with Wat(0,x, k) = A(x, k), swe/t i/tai Wat G C 1 ([0, +oo), L 2 (M m )). Wat is owen 6y t/ie series 
expansion 

W N (t, x, k) =e" Sjvi A(x + ik, k) 

+ Yl I ds(n) I d p (n] ft ^) A ( x + tk + E w» k + E Pi) • (34) 

n >l JD n{t) J j = l j = l j-l 

Here, D n (t) = {s^ = (si, . . . ,s n ) : < s n < • • • < si < i}, p (n) = (pi, . . . ,p„), and 

Sat = / dpo"(p)- 

«/|p|>l/W 

Let W be a weak solution of the radiative transfer equation (|18|) with a null initial condition. 
Using the test function Wn given by Proposition 16.31 let us consider Wjv(^x, k) = Wn(T — 
t,x, k), for T > 0. Then, Wn satisfies 

d t W N (t, x, k) + k • V x W N {t, x, k) + C N W N (t, x, k) = 

with Wjsr(T, x, k) = A(x, k). Consequently, we have following lemma. 
Lemma 6.10 Vt G (0, T), we /iawe 

^<^(t), ^(t)) L2(R2d) = WaK*) + k • V x WV(i) + CW N (t)) L2(R2d) 

= (W(t),(C-C N )W N (t)) L2iR2dr 

As a result, 

(W(T),X) L2 2d =(W(T),W N (T)) L2 2d = f ds(W(s),(C- C N )W N (s)) L2 2d 

Finally, it remains to show the following lemma. 
Lemma 6.11 

sup ||(£ - C N )W N {t)\\i>Q*«) <Cx[ 
te[o,T] J\p\<i/n 0(p) 

This last lemma gives us the weak uniqueness. In fact, we have / <ip-Ro(p)|p|/fl(p) < 
and then, 

lim f d P ^M=0. 

Consequently, VA G C^(R 2d ) and VT > 0, (W(T), A) i2(B;2d) = 0. 

Proof (of Propostion 16. 3|) Let us rewrite (|33p using the finite scattering coefficient Eat- 
Then, We are looking for a solution of 

dtW N (t,x,k) =k-V x W J v(t,x,k)-E JV W i v(t,x,k)+ /" d P (j(p)Ty 7 v(t,x,k + p) 
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This equation can be recast by integrating as follows 

W N {t,x,k) = e- Sjv *A(x + tk,k) + [ ds f c2 P (j(p)e~ ( *~ s)s, Wjv(s,x + (i-s)k,k + p), 

JO J\p\>l/N 

and then the series expansion (|34|) is obtained by induction. Let us write 

W N (t,x,k) = £Wft(t,x,k) 



n>0 



where are the terms of the serie expansion (|34p . Using changes of variables we show that 
Vi > 0, W N {t) G L 2 (M d ) and 

X) BU PlTO( t )Hi a (tt M ) < +0 °- 



n>0 * 



Then W N G C°([0, +00), Z/(R 2d )). Moreover, 

^ sup ||k • V x W#(i)|| L 2( R 2d) < ||k • V x A|| i2 ( R 2d) + CTS at || V X A || L 2 (R 2d) < +00, (35) 



n>0°<*< T 



since Rq is assumed to be with compact support, and 

SU P II / dpa(p)W%(t, ■, ■ + p)|| L 2 (R2d) < E N \\X\\ L2(MM) < +00. (36) 

n>0 * 



Proof (of Lemma I6.10|) First of all, let us note that equation (|18p is valid for all A G 

C(f(M 2d ). However, since sup t > HW^H^d) < +00 and W G C°([0, +00), (£ r , d Br )), we 
also have 



(W(t), A(t)) L2(R2d) = (W(0), A(0)) L2(R2d) + / ds(W(s), d t X(s) + k • V X A( S ) + CX(s)) L2 {R2dy 



VA G C 1 ([0,T],C x {R 2d )), for which sup t ||A(t)|| £ 2 (R 2d) + sup t ||a t A(t)|| L 2 (R 2d ) < +00. Conse- 
quently, we have 

(w(t), wm) L 2 (R2d) =(w(o), W%(0)) L2(R2d) 

+ jT ds(W(s),d t WMs) + k • V x W#(s) + CW$(s)) La 
where W%(t) =W%{T-t). First, let us note that 

W#(t,x,k) = / ds [ dpa(p)e-^ Sjv ^- 1 (s,x+(t-s)k,k + p), 

J J|p|>l/A r 

so that, thanks to ([33]) and ([55]) . we obtain 



S „! U F m ll S t^(*)IU 2 (M^) + I] „! U P m H k • V x^(0llL2(R2d) < +00- 

Moreover, 



n>0°<*< T n>0°<*< T 



CWjtft, x, k) = y dpa(p) ^ dnp • V k W#(i, x, k + up) 
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and then 

]T sup \\£Wm\\mR^)< I rf P ^ P) ' Pl (V||V X A|| L2(RM) + ||V k A|| L2(R2d) ) < +oo. 

Consequently, we obtain the desire result 

(W(t), W N (t)) L2(R2d) =(W(0), W N (0)) L2(R2d) 

+ J* ds(W(s),d t W N (s) + k • V x W N (s) + CW N (t)(s)) l 
=(W(0),W N (0)) L2{M2d) 

+ f ds{W{s),{C-C N )W N {s)) L2 



L 2 



□ 

Proof (of Lemma 16. lip We have 



(£ - £ N )W N (t,x,k) = [ dpa(p) [ dup ■ VkW N (t,x,k + up) 

J\p\<i/N Jo 



and then 



sup ||(£ - C N )W N (t)\\ L2{R2d) < I dp ^o(p)|p| / T||VxA|| + || VkA || 

)<t<T ^ ; J\v\<HN P V V ' 



0<t<T J\p\<l/N Q{P) 

□ 

7 Proof of Theorem 13.11 



First, we let us compute the solution of (|T8|) . 

Lemma 7.1 VW G L 2 (R 2d ), the unique weak solution uniformly bounded in L 2 (M. 2d ) of (fTSj) 
is given by 

W(t, x, k) = ^-Jp J dydqe^^ exp ( J* du^(q + uy)) W (y, q + ty). (37) 
where ^ defined by (|21|). 

Proof (of Lemma El]) First, let W G Cg°(lR 2d ) and let 

W(t, y, q) = exp ( J du^{q + uy)) Wb(y, q + ty). 

W is the inverse Fourier transform of W. Then, 

d t W = y ■ V q VF + -V(q)W, 

and taking the inverse Fourier transform of this last equation we obtain that W is a classical 
solution, and then a weak solution of (|18|) uniformly bounded in L 2 (M. 2d ). In fact, 

\\W(t)\\v mM) < \\W \\ LHR 2 d) (38) 

since 

|exp ^ J du^(q + uy)) = exp ^ — J du J dpcr(p)(l — cos(p • (q + uy)))) < 1. 

Finally, using (|3"8"j) . the uniqueness of (fTHj) . and by density, we obtain that ([3T|) is the unique 
solution of (|18|) . That concludes the proof of Lemma 17.11 □ 
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Now, let us show the key lemma of the proof of Theorem 13.11 This lemma is a consequence 
of assumption (fT3"j) . 

Lemma 7.2 Let \i > and t £ (3/(2//), 2///). There exists M( )(U > and C^e^ > 0, such that 
V|q| > and V|y| > 

exp(^ d«*(q + ^y)) < exp(-C tjeiM |q| e )l ( | q |> At4 | y | ) + exp(-C tiejA1 |y| e )l ( | q | <M | y | ) . 

Proof (of Lemma 17. 2|) This lemma is a consequence of assumption (|13p . First, we have 

/ du^(q + ufiy) = - / du dpa(p)(l - cos(p • (q + ufiy))) 
Jo Jo J 

-~ dU I | P |<l/(2|q|) d P fJ (P)( 1 - COS (P • (q + *W))) 
J ° J |p|<l/(2^|y|) 

-~i dU l IPl<V(2|q|) ^(P)!? • (*1 + U ^)| 2 



Then, 

rt t 3 
/ du|p • (q + ufiy)\ 2 = t\p ■ q| 2 + fit 2 p ■ q p • y + /i 2 — |p • y[ 

JO o 



2 



= (* - m^-)Ip • q| 2 + {v 2 j - a^Ip ■ y| 2 + m^-Ip • (q + y)l 2 
> t(i - ^)|p • q| 2 + /^ 2 (m~ - j)Ip • yl 2 ' 



>0 >0 

since t £ (3/(2/i), 2//i). Consequently, 



f du^in + liyuy) < -C t ^{\p ■ q| 2 + |p • y| 2 ), 
Jo 



and 



exp ( I du*(q + u^yf) < exp(-Cf iM j | p | <1/(2 | q |) dp<r(p)|p • ( l\ 2 ) 1 (H>^\y\) 

° ' | P |<i/(2M|y|) 

+ eX P (-C t)/ , y |p|<l /(2 |q|) ^(p)|p-y| 2 )l(|q|<M%|)- 



|p|<l/(2|tt|y|) 

Moreover, according to (fTBj) . there exists Mt i(U large enough so that V|q| > M t ^ and V|y| > 



expl / du^{q + ufiy)) < exp I - C t ^\q\ u I dp \p • -M" j l(| q |>^| y |) 



/o 



/ dp 

J|p|<l/2 


1 


P • 






\p\d+e 


M 


i 2 ) 


/ ^Pl 

/|p|<l/2 


l 


P 


— 1 




p ci+6» 


|y| 


2 > 



(39) 



Finally, thanks to a unitary transformation on § 



d-i 



d VT-uTa |p ■ rrl 2 = / l u • v| 2 dS(u) 



|p|<l/2 ~ |p| d+e |Z| 1 7§rf-l' 

where v € S d_1 is fixed. □ 
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The following lemma shows the transport of regularity from the variable k to the variable 
x through the transport term in (|18|) . This transport of regularity is a direct consequence of 
Lemma 17.21 

Lemma 7.3 Let W be the unique weak solution uniformly bounded in L 2 (M. 2d ) of (|18p . We 
have Vii > 

W G J L°°([ii,+oo), p| H k (R 2d )) 

Proof (of Lemma [773]) Let h > 0, t > h, and rj > 0, such that t £ (3/(2//), 2///). Then, 
using a change of variable, we have 

\\W(t)\\ 2 Hk{R2d) <C, j q| > Mtft dydq(|q|* + |y| fc )e/o ***(«»+w)|Wo(y,q + ^*y)| 2 

\y\>Mt' tfl 

+ C » U>M t „ dyd^elo ^+^y)\W (y, q + Vty)\ 2 
\y\<Mt^ 

+ C m / q |<M t „ dydq|y| fc e/o rf «*(^)|^ ( y , q + #y)\* 

\y\>M t ',ij, 

+ C m Li 

\y\<Mt',p 

The first and last terms are finite since Wq E L 2 (R 2d ) and thanks to Lemma 17.21 The second 
and third terms are finite using (|39|) depending on whether q or y is bounded. □ 

Finally, from Lemma 17.31 we have the continuity of W. 

Lemma 7.4 Let W be the unique weak solution uniformly bounded in L 2 (R 2d ) of (|18|). We 

have 

W G C°((0,+oo), Pi H k (R 2d )). 

k>0 

Proof (of Lemma 17.4"]) Let t± > and (t n ) n be a sequence which converges to t\, Then, 
there exists N and to > such that Vn > N, t\ At n > to, and 

||W(t n )-^i)ll^ (R 2 d) <C fc | dy(iq(|y| fc + |q| fe )|^(to,y,q + tiy)l 2 

x |i_X* (q+(u "* o)y) | 2 

+ C k J dydq(|y| fe + |q| fe )|W(i ,y 5 q + tny)- W(*o, y, q + *iy)| 2 - 
Then, we can conclude from Lemma 17.31 and the dominated convergence theorem. □ 

8 Proof of Proposition 14.11 

First, let Wo £ Cq°(1R m ) and let (Lt)t be a Levy process with characteristic exponent VP defined 
by (|2ip and infinitesimal generator given by C defined by (|19|) . and so that (-Lt)t is a Markov 
process on Co i0 o(K 2d )) the space of continuous functions which tend to at infinity. Moreover, 
let us note that C^{R 2d ) C £>(£). 

According to assumption (fT3"]) . the fact that a has a support included in a compact set, and 
[351 Lemma 31.7], the process defined by 

X t = (x-t\n- j\„da,k + Lt) (40) 
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is a Markov process on Co i0 o(I^ 2<i ) with infinitesimal generator given by 

L(p(x, k) = — k • VxV( x ; k) + Cip(x, k). 
As the result, The function defined by 

W(t, x, k) = E [Wo (x - tk - J L s ds, k + L 
is a classical solution of 

d t W + k • VW = CW. 
Moreover, if pit, dx.dk) is the transition function of the Markov process (/q L s ds, L 

\\W(t)\\ 2 L2(R2d) = J dxdk\ J Ty (x-tk-q,k + q)p(t,dy(iq)| 2 

< J p{t, dydq) J dkJdx\W (x- tk- q,k + q 

< IIWo|i! 2(KM) , 



tit, 



(41) 



(42) 



so that W is uniformly bounded in L 2 (M. 2d ) and is the unique weak solution of (|18|) , However, 
(|4"2"j) shows that (jUJ) can be extend to W G L 2 (M 2d ). Then by density, (jH]) is the unique 
weak solution uniformly bounded of (fT8"j) with Wo G L 2 (M 2c( ). That concludes the proof of 
Proposition 14.11 



9 Proof of Theorem 15.1 

First, Let us recall that 

supsup WWWWv^ < \\W \\ L 2 {W 2 d) , (43) 

V t>0 

where W v is the unique solution uniformly bounded in L 2 (M. 2d ) of the radiative transfer equation 
(fT8"j) associated to ^(p). 

The proof of Theorem 15,11 uses the probabilistic representation obtained in Proposition 14.11 
The following proof is in two step. First, we prove the relative compactness of (W T, )„, afterward 
we identify all the subsequence limits. The proof of the relative compactness of (W v ) v is based 
on the probabilistic representation obtained in Proposition 14. 1\ but where is associated 

to the jump measure cr^p^p. 

Lemma 9.1 V(x,k) € R 2d , (X^) v , defined by (@QD, is tight is on V([0, +oo), R 2d ). 

Proof (of Lemma 19. ip To prove this lemma we will use the Aldous tightness criteria [11]. 
Let us begin with the tightness of {L^) t . First of all, let us decompose (L^) t according to the 
size of its jumps, 

L\ = Ll>" + L 2 " = [ pN^(t,dp)+ [ pN*(t,dp), 
J\p\<i J]p|>i 

where dp) is a random Poisson measure with intensity measure ^(p)^. Let us note 

that (Ll' v ) t is a Levy process with compactly supported jump measure o"' 7 (p)l|p|< 1 (ip, so that 

sup sup EflL^I 2 ] < +oo, 

V 0<t<T 
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thanks to ([25]) and [U Theorem 2.3.7]. Moreover, (Lf' 1] )t is a compound Poisson process [H 
Theorem 2.3.9]. Then, there exist a poisson process (Nt)t with intensity 1 and a sequence 
(Z%) n of iid random variables with distribution cr 11 (jp)l\ p \>idp, such that 

N(t) 
j=0 

in distribution with Zq = 0. Moreover, (Nt)t and (Z%) n are independent. First, 



( sup |Z,J?| > M) < P( sup |L^ ,r? | > Af/2) +P( sup > Af/2). 

0<i<T 0<t<T 0<t<T 



Since 



Lj'"-t/ P^(p)dp 
«/|pI<i 

is a martingale, and for Af sufficiently large, 

P( sup \L x t *\ > Af/2) < P( sup - t [ pa r '{p)dp > Af/4 

0<t<T V 0<i<T f|p|<l 

; S u P E[|L^| 2 ] + cr(/ — ^L_rfp 
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according to the Doob's maximal inequality. Moreover, 



P( sup \L;> V \ > Af/2) < V P( V |Z?| > Af/2 JV(T) = n)P(JV(T) 

()</' 7 ' •' 



n>0 j=0 
n 



Then, 



< ^^P(|2j| > Af/(2n))P(iV(T) = n) 

n>0j=l 

< V n [ a v (p)dp¥(N{T) = n) 

«'[p|>Af/(a»») 

< ^W^]. 
supP( sup > Af) = 0(—7r) as Af -)• +oo. 



i) 0<t<T ^Af 61 ' 

Consequently, VAf ' > 

ft 



( sup / L^ds > Af ) < P( sup |2#| > Af') +P( / m\ds > Af, sup |L?| < Af'). 

V 0<t<T JO J V 0<t<T ' K h 0<t<T ' 



and then, 



lim supP( sup \X^(t)\ > Af) = 0. 

A/^+oo r; 0<t<T 



Let r be a discrete stopping time relatively to the filtration of (X 11 )^, and bounded by T. 
According to the stationarity of (L \ )t, V/i > and W > 0, we have 



T+fJ, 



P(|X"(r + fj.) - X"(r)| > u) < P(|L"(r + /x) - L"(r)| > u/2) + F^j \Ul\ds > v/2) 

< P(|L"(/i)| > v/2)+F(J^m\ds > v/2). 
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First, it is clear that MM' > 

P( r\Ul\ds > v/2) < p( sup \m > M') + p( r\Vl\ds > i//2, sup \L%\ < M 1 
Wo ' v o<t<r ' Wo o<t<T 

Moreover, 

P(|£"(M)| > f/2) < ]P(|^ 1,7? (m)| > + P(|L 2 ' r '(^)| > i//4) 

However, E[N(fi) 1+e ] < E[A%) 2 ] = \i + ^i 2 , so that 

lim supsupPQX^r + n) - X t '(t)\ > v) = 0. 

That concludes the proof of Lemma 19.11 □ 

Lemma [9 . 1 1 implies the relative compactness of (W ri ) ri . Now, let us identify all the accumu- 
lation points of (W v ) v . First, let us denote 

C~<p(k)=Jdpo™(p)( ( pQz + p)-<p(k)) and ^(p) = _^-_. 

We have the following lemma which shows the convergence of the transfer operator C^. 
Lemma 9.2 We have VA G Cg° 



lim||£n-£°°A|| L2(Rd) =0. 
Proof (of Lemma 19. 2D This lemma is a consequence of 

^£V\ _ £ 0O A|| L 2( R 2d) < C(||A|| L 2( R 2d) + ||VkA||) 

x([ dp|p||a"(p)-<7°°(p)| + / dp\a«(p)-a°°(p)\), 

and the dominated convergence theorem, because of (|24|) and (|25p . □ 
Now, let us consider r\ n a sequence, which goes to as n — > +oo, and such that V(t, x, k) 

lim W^(t,x,k) = W°°(t,x,k). 

n— >+oo 

Lemma 9.3 VF°°(i,x, k) is a weak solution uniformly bounded in L 2 (M. 2d ) of the radiative 
transfer equation 

dtW 00 + k • V X W°° = C^W 00 with W°°{0, x, k) = W (x, k). (44) 



Proof (of Lemma 19. 3j) First, by Fatou's Lemma, we have Vi > 

\\W°°(t)\\L*<X") <Um n inf||W»(t)|| i2(R2d) < ||Wo|| i2( 
From this result, using a molifier, and the dominated convergence theorem, we have that 

Vt>0, lim W Vn (t) = W°°(t) weakly in L 2 (R 2d ). 

n— >+oc 

According to (|43p . Lemma 19.2} the dominated convergence theorem, and taking the limit in 

(W*(t), A) L2(R2d) - (Wo, A) L2(R2d) - f (W^(s),k ■ V X A + C^\) L2 ,^ds = 
VA G Cg°(M 2d ), we get the desired result. □ 
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The following lemma concludes the proof of Theorem 15.11 

Lemma 9.4 Equation ([34"]) admits a unique weak solution uniformly bounded in L 2 (M. 2d ), which 
is given by ([27|) . 

Proof (of Lemma 19. 4|) First of all, let us recall the following Levy-Khintchine formula [35|. 
Thorem 14.14] 

/ - 1} = - |q| ^ r(1 - e) jL dS[u)lei ' u|9 ' 

which give VA G C^(M. 2d ) 

£°°A = -a(9)(-A) e / 2 X. 
To show the uniqueness of Equation ([33]) , let us defined 

A(t, x, k) = ^ J dydqe i ^y+ k ^ exp ( - a(6) jT^q + uy\ e ) A (y, q + (T - t)y), 

with A G Cg°(M 2d ), and let TV be a solution of (03]) uniformly bounded in L 2 (R M ) with 
W(0) = 0. Then, since V(C°°) = H e (M. d ), we have 

(W(T), X ) L2(R2d) = (W(T), A(T)) L2(R2d) = f T ds(W{s) : d t X + k • V X A + £°°A) i2(R2d) = 0. 

Moreover, it is clear that (]27|) is a weak solution for Wo G Co°(]R 2rf ). Consequently, by density 
of Cg°(R 2(i ) in L 2 (R 2d ), ([27]) is also a weak solution uniformly bounded in L 2 (R 2d ) with W G 
L 2 (R 2d ). □ 

Conclusion 

In this paper we have analyzed the asymptotic phase space energy density of a solution of the 
random Schrodinger equation with long-range correlation properties. The phase of a solution 
of the random Schrodinger equation in the same context has been studied in [TU], and the 
phase and the phase space energy density for the random Schrodinger equation with rapidly 
decaying correlations has been studied in [H] E] [TU] . In the case of rapidly decaying correlations 
the phase and the phase space energy density evolve on the same scale, and the asymptotic 
evolution of the phase space energy density is given by a radiative transfer equation. In the 
other hand, for random perturbations with slowly decaying correlations, the phase and the 
phase space energy density do not evolve on the same scale anymore. However, the scale 
of evolution of the phase space energy density in the context of rapidly and slowly decaying 
correlations are the same, and their evolutions are given by radiative transfer equations of the 
same form, but nevertheless with thorough differences. In the context of perturbations with 
long-range correlations, a scattering coefficient cannot be defined because of a nonintegrable 
singularity in the radiative transfer equation. However, this singularity is responsible to a 
regularizing effect on this radiative transfer equation which cannot be observed in the case 
of rapidly decaying correlations. Finally, we have derived an approximation of the radiative 
transfer equation involving a fractional Laplacian, which gives us a simpler model of phase space 
energy distribution which corresponds to the long space and time diffusion approximation of 
the radiative transfer equation. Moreover, this model permits to exhibit a decay of the phase 
space energy distribution which obey to a power law with exponent lying in (0, 1). 
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